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A GRAPHICAL  METHOD  OF  EVALUATING 


CERTAIN  CRYSTALLOGRAPHIC  STRUCTURE  FACTORS 

by 

H . .T . Grenville-W ells 

Laboratory  for  Insulation  Research 
Massachusetts  Institute  of  Technology 
Cembridge,  Massachusetts 

Abstract:  A graphical  method  is  described  for  evaluating  structure  factors  of 

n 

the  form  F(hk)  = £ * (hk)  cos  2ir(hx  + ky  ) and  electron  density  expres- 

o 

sions  of  the  form  p(X,  Y)  = £ £ F(hk)  cos  2ir(hX  + kY).  The  method  de- 

r H K 

pends  upon  the  fact  that  if,  in  a [ooij  projection  of  a unit  cell  the  normal 
to  the  planes  (hkO)  is  graduated  with  the  appropriate  cosine  function  giving 
the  phases  of  points  along  the  normal;  then  a circle,  with  its  center  at 
(x/2,  y/2),  passing  through  the  origin  will  cut  the  normal  to  F(hkO)  at  the 
value  of  cos  2ir(hx  + ky). 


Introduction 

The  construction  discussed  below  provides  a rapid  graphical  method  for 
evaluating  Fourier  summations  involving  trigonometric  functions  of  the  form 
cos  2ir(hX  + kY),  or  functions  such  as  cos  2irhx  • cos  2-rrky  which  can  be  reduced 
to  this  form. 

It  is  intended  to  apply  to  two-dimensional  projections  of  a crystal  struc- 
ture, and  the  nomenclature  used  below  reiers  to  a [00l]  projection.  It  can 
be  used  either  on  the  true  projection  of  the  unit  cell,  or  with  the  atomic  co- 
ordinates referred  to  a square  projection,  a procedure  which  has  some  ad- 
vantages to  be  discussed  in  Section  3.  The  real  lattice  and  not  the  reciprocal 
lattice  is  used  in  these  particular  applications. 
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1 . Basic  Construction  for  Evaluating  Structure  Factors 
Consider  the  [001]  projection  of  the  unit  cell  (Fig.  1).  AB  is  the  plane 
(hkO)  so  that  OA  = a/h,  OB  = b/k.  OP  is  the  plane  through  O parallel  to  AB. 
The  contribution  to  the  structure  factor  F(hk)  from  the  atom  at  (x.  y)  can  be  ob- 
tained  from  a Bragg-Lipson  chart,  and  the  line  PQ  can  be  graduated  from 
such  a chart.  The  particular  value  of  this  construction  lies  in  the  fact  that 
since  PQ  must  be  perpendicular  to  OP,  points  such  as  P will  lie  on  the  circle 
OPQ  for  all  planes  (hkO).  This  is  the  circle  with  its  center  at  (x/2,  y/2)  which 
passes  through  the  origin.  In  Fig.  2,  ON  is  the  normal  through  the  origin  to  the 
planes  (hkO),  and  it  can  be  seen  that  OR  is  equal  and  parallel  to  PQ.  Hence,  if 
ON  instead  of  PQ  were  graduated  from  the  appropriate  Bragg-Lipson  chart,  the 
value  of  f ^ ( hk)  cos  2ir(hxr  + kyf)  would  be  the  reading  at  the  point  where  the 
circle  with  center  (xr/2,  y^/2)  intersects  ON;  ON  will  be  called  the  "F(hk) 
line". 

Generalizing  the  construction,  it  can  be  seen  that  for  any  (hkO)  plane,  the 

structure  factor  F(hkO)  will  be  the  sum  of  the  readings  on  the  corresponding 

F(hk)  line  where  it  is  intersected  by  the  n circles  corresponding  to  the  n atoms 

n 

in  the  projection,  i.e.,  F(hk)  = £ f (hk)  cos  2n-(hxr  + kyf).  The  method  is  il- 

o 

lustrated  in  Fig.  3. 

2.  Basic  Construction  for  Evaluating  Electron  Density 
By  analogy  with  the  preceding  construction,  if  a set  of  normals  to  the 
planes  (hkO)  are  graduated  as  before  from  a Bragg-Lipson  chart,  but  instead 
of  carrying  the  values  of  the  atomic  scattering  factor  fQ(hkO),  they  are  given 
the  observed  values  of  F(hkO)  (Fig.  4),  a circle  whose  center  is  at(X/ <i,  \/ l.) 
will  cut  the  F(hjk1)  line  at  A,  which  will  give  the  value  of  Ffhjkj)  cos  2tt 


1)  W.L.  Bragg  and  H.Lipson,  Z.  Krist.  A95,  323  (1936). 
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Fig.  1.  Basic  construction  for  structure  Fig.  2.  Use  of  normal  to  (hkO)  planes 


factors  of  the  form  F(hk)  = £f  (hk)  cos 

o 

2ir(hxr  + kyr). 


to  represent  ; (hkO). 


(hjX  + kjY)  and  so  on.  Hence,  the  sum  of  the  readings  on  the  circle  OABCD, 

which  will  be  called  the  "^>(X,  Y)  circle,"  will  be  F(hjk, ) cos  2Tr(hjX  + k^  Y)  + 

F(h  k ) cos  2tr(h  X + k Y)  + for  all  the  observed  values  of 

2 2 2 2 H K 

F(hkO),  i.e.,  the  sum  will  be  2 2 F(hk)  cos  2ir(hX  + kY),  which  is  the  value 

O O 

of  the  electron  density  p[X,  Y)  at  the  point  (X,  Y). 

3.  Use  of  a Square  Projection 


There  are  two  main  advantages  in  the  use  of  a square  projection:  (1)  The 
directions  and  graduations  of  the  F(hk)  lines  will  then  be  the  same  for  all  cry- 
stals, so  that  a master  set  can  be  prepared  for  all  projections  instead  of  a new 
set  being  necessary  for  each  one  (except  in  the  case  of  axial  projections  of  cubic 
crystals  or  [001]  projections  of  tetragonal  crystals,  which  are  already  square), 
and  only  the  numerical  values  of  f(hk)  or  F(hk)  will  have  to  be  inserted  for  each 
projection.  (Z)  li  Bragg-Lipson  charts  are  not  available,  the  F(hk)  lines  can  be 
graduated  with  very  little  labor  by  making  use  of  the  construction  shown  in  Fig.  5 
if  a square  projection  is  used. 


-5- 


T.R.65 


If  a circle  drawn  with  its  center  at  the  point  ( 1 / 2 , l/Z)  intersects  the  F(hk) 
line  at  P,  OP  has  to  be  divided  into  (h  + k)  parts;  i.e.,  in  Fig.  5,  F(4,  2)  is  di- 
vided in  4 + 2 = 6 parts.  If  the  charts  are  drawn  on  tracing  paper,  a piece  of 
graph  paper  can  be  conveniently  used  as  a proportional  divider  to  graduate  the 
F(hk)  lines.  A disadvantage  is  that  the  crystal  projection  is  distorted. 


sily  summations  of  the  form  P!X,  Y)  = 
H K 

£ £ F(hk)  cos  2ir(hX  + kY). 

o o 


Fig.  5.  Graduation  of  F(hk)  lines  using 
a square  projection. 


4.  Extension  to  Structure  Factors  of  the  Form 
n 

F(hk)  = £ f (hk)  cos  2iThx  • cos  2irky 
o 


The  expression  fr(hk)  cos  2-irhx^  • cos  2-sTkyr  can  be  rewritten  in  the  form 

1/2  £t'r(hk)  cos  2ir(hxr  -t-  kyf)  + fr(hk)  cos  2Ti(hxr- kyr)  J . ...  (1) 

This  is  equivalent  to  placing  half  an  atom  at  the  point  (xr,yr)  and  l/ 2 an  atom  at 
(x  ,y  ).  Then  in  Fig.  7 the  value  of  Eq.^l)  will  be  one  half  the  sum  of  readings 
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at  A and  B.  Although  an  atom  at  (xr>  y ) is  actually  located  at  A (Fig. 6)  so  that 
its  contribution  to  F(hk)  occurs  at  P,  it  can  be  seen  that  if  O Q is  drawn  paral- 
lel to  OP,  in  which  case  OQ  = OP,  the  reading  at  Q will  be  equivalent  to  the  read- 
ing  at  R,  and  therefore  the  atom  at  (xr«  yr)  can  be  placed  at  B and  the  construc- 
tion confined  to  one  quadrant.  Atoms  with  co-ordinates  (x,  o)  will  be  unaffected  and 


Fig.  6.  Equivalence  of  atomic  positions  (x,  y)  and  (x,  1 — y). 


7.  Construction  used  for  structure  factors  of  the  form  F(hk)  = ^ fr(hk.) 


cos  2irhx  • cos  2irky  . 
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5 . Cases  Where  hkO  Reflections  Are  Required  in  Addition  to  hkO  Reflections 
The  construction  as  discussed  so  far  has  been  confined  to  a single  qua- 
drant, i.e.,  to  F(hk)  lines  where  h and  k are  both  positive.  For  some  projec- 
tions F(hE)  reflections  will  not  be  equivalent  to  F(hk)  reflections,  and  will  have 
to  be  included  separately  in  the  summations.  It  can  be  seen  in  Fig.  8 that  the 
circles  with  centers  at  (0,  1 / Z) , ( l/ Z , 0)  and  ( 1 / Z , 1 / Z)  define  a region  within 
which  all  other  (x,  y)  circles  will  lie. 

This  region  does  not  extend  into  quadrant  (3),  and  hence  extension  of  the 
F(hk)  line  will  not  be  necessary;  but  it  does  fall  partly  in  quadrants  (Z)  and  (4), 
so  that  the  F(hk)  line  will  have  to  extend  through  the  origin  into  quadrant  (4). 


6.  Circles  Tangential  to  F(hk)  Lines 
An  (x,  y)  circle  will  be  tangential  to  an  F(hk)  line  only  if  (x,  y)  lies  on  a 
line  through  the  origin  perpendicular  to  the  F(hk)  line,  i.e.,  if  (x,  y)  lies  in 
the  (hk)  plane.  Hence  it  will  make  a full  ^.(hk)  contribution  to  F(hk). 
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7 . Practical  Considerations 

There  are  several  limitations  on  the  practical  usefulness  of  the  construc- 
tion. 

(a)  If  the  F(hk)  lines  are  graduated  with  values  of  the  atomic  scattering 
factor  f(hk),  this  becomes  confusing  when  several  kinds  of  atoms  are  present. 

For  two  types  of  atom  two  independent  sets  of  numbers  can  be  used  (one  on 
each  side  of  the  F(hk)  line),  but  more  would  be  impracticable.  It  would  be  pos- 
sible to  prepare  a separate  chart  for  each  species  of  atom,  or,  in  cases  where 
only  F(hk)  with  h and  k positive  occurs,  i.e.,  in  which  only  one  quadrant  is 
used,  to  construct  a multiple  chart  (hypothetical  example  given  in  Fig.  9)  with 

different  f(hk)  values  in  each 
quadrant.  The  actual  number 
of  graduations  used  on  the  F(hk) 
line  will  depend  largely  on  the 
scale  chosen  for  the  figure. 

If,  for  example,  the  sym- 
metry of  the  projection  is  such 
that  only  atoms  in  the  section 
bounded  by  a/4  and  b/4  need  be 
considered,  then  the  use  of  a 
25  cm  square  jvould  mean  that 
I niiii  - 0 . 001  of  the  cell  edge. 

In  this  case  the  largest  interplanar 
spacing,  d,^  will  be  100  cm  on 

the  chart,  i.e.,  360°  = 100  cm,  and  values  of  f(hk)  cos  2ir(hx  + ky)  could  be  in- 
serted at  1°  intervals.  This,  however,  is  an  extreme  case,  and,  in  the  examples 
jq  fn »•  constructed,  values  of  f(hk)  cos  ?i?(bx  v kv)  have  been  inserted  at  15  inter  — 


Fig.  9.  Composite  chart  for  structure  con- 
taining several  types  of  atom. 
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vals  from  Table  1,  showing  these  values  for  f(hk)  = 1,  2 . . . . 99;  intermedi- 
ate values  have  been  interpolated. 

(b)  If  all  observed  F(hk)  lines  are  inserted  in  the  diagram,  it  becomes 
hopelessly  crowded.  It  would  be  possible  to  indicate  the  positive  and  negative 
sections  of  a large  number  of  F(hk)  lines,  but  use  of  numerical  values  would 
be  impracticable. 

(c)  It  is  very  difficult  to  show  several  orders  of  F(hk)  on  the  same  line 
clearly.  These  are  serious  handicaps  as  regards  use  of  the  construction  for 
electron  density  summations,  when  contributions  from  all  F(hkO)  values  must 

be  included,  but  can  be  circumvented  in  the  case  of  structure-factor  calculations, 
and  even  turned  to  advantage. 

8.  Introduction  of  the  Temperature  Factor 

Consider,  for  example,  a crystal  in  which  two  kinds  of  atoms  A and  B 

are  present,  and  suppose  that  the  atomic  scattering  factors  of  these  atoms 

at  rest  can  be  represented  by  the  curves  A and  B in  Fig.  10.  Suppose  that  C 

and  D,  respectively,  represent  the  smallest  and  largest  observed  values  of 

(sin  8/X.).  The  range  CD  is  now  divided 

into  sections  like  CEFG  so  that  f (C) 

o 

and  f (G)  do  not  differ  from  f (M)  by 
o o 

more  than  5 percent. 

A chart  such  as  described  in  Section 
1 can  now  be  constructed  with  only  those 
F(hk)  lines  for  which  l/2d  lies  in  the 
range  EF,  say,  8 F(hk)  lines,  and  this 
will  be  quite  manageable.  The  fQ(A)  and 
f (B)  values  are  used.  Whatever  the 
value  of  the  temperature  factor,  the  ratio 


it;  „ in  a * 

M.  » A W . .A.Wl.l 

curves. 


f o rf  r\  1 
*6  


I 
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Table 

1. 

Values  of  A cos  2ttR 

, for 

A = 1 

• 

• • 

99;  R 

ii 

o 

o 

X 

1/24 

-►  1/4 

• 

0° 

15° 

30° 

45° 

60° 

75° 

90° 

0° 

, - G 
1 0 

OJ 

o 

0 

45° 

60° 

75° 

90° 

1 

1 

1 

1 

0 

0 

0 

51 

49 

44 

36 

26 

13 

0 

*> 

2 

2 

1 

1 

1 

0 

52 

50 

45 

37 

26 

13 

0 

3 

3 

3 

2 

2 

1 

0 

53 

51 

46 

37 

26 

14 

0 

4 

4 

3 

3 

2 

1 

0 

54 

52 

47 

38 

27 

14 

0 

5 

5 

4 

4 

2 

1 

0 

55 

53 

48 

39 

28 

14 

0 

6 

6 

5 

4 

3 

2 

0 

56 

54 

48 

40 

28 

14 

0 

7 

7 

6 

5 

4 

2 

0 

57 

55 

49 

40 

28 

15 

0 

8 

8 

7 

6 

4 

2 

0 

58 

56 

50 

41 

29 

15 

0 

9 

9 

8 

6 

4 

2 

0 

59 

57 

51 

42 

30 

15 

0 

10 

10 

9 

7 

5 

3 

0 

60 

58 

52 

42 

30 

16 

0 

11 

11 

10 

8 

6 

3 

0 

61 

50 

53 

43 

30 

16 

0 

12 

12 

10 

8 

6 

3 

0 

62 

60 

54 

44 

31 

16 

0 

13 

13 

11 

9 

6 

3 

0 

63 

61 

55 

45 

32 

16 

0 

14 

14 

12 

10 

7 

4 

0 

64 

62 

55 

45 

32 

17 

0 

15 

14 

13 

11 

8 

4 

0 

65 

63 

56 

46 

32 

17 

0 

16 

15 

14 

11 

8 

4 

0 

66 

64 

57 

47 

33 

17 

0 

17 

16 

15 

12 

8 

4 

0 

67 

65 

58 

47 

34 

17 

0 

18 

17 

16 

13 

9 

5 

0 

68 

66 

59 

48 

34 

18 

0 

19 

18 

16 

13 

10 

5 

0 

69 

67 

f'0 

49 

34 

18 

0 

20 

19 

17 

14 

10 

5 

0 

70 

68 

61 

49 

35 

18 

0 

21 

20 

18 

15 

10 

5 

0 

71 

69 

61 

50 

36 

18 

0 

22 

21 

19 

16 

1 1 

6 

0 

72 

70 

62 

51 

36 

19 

0 

23 

22 

20 

16 

12 

6 

0 

73 

71 

63 

52 

36 

19 

0 

24 

23 

21 

17 

12 

6 

0 

74 

71 

64 

52 

37 

19 

0 

25 

24 

22 

18 

12 

6 

0 

75 

72 

65 

53 

38 

19 

0 

26 

25 

23 

18 

13 

7 

0 

76 

73 

66 

54 

38 

20 

0 

27 

26 

23 

19 

14 

7 
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77 

74 

67 

54 

38 

20 

0 

28 

27 

24 

20 

14 

7 

0 

78 

75 

68 

55 

39 

20 

0 

29 

28 

25 

21 

14 

8 

0 

79 

76 

68 

56 

40 

20 

0 

30 

29 

26 

21 

15 

8 

0 

80 

77 

69 

57 

40 

21 

0 

31 

30 

27 

22 

16 

8 

0 

81 

78 

70 

57 

40 

21 

0 

32 

31 

28 

23 

16 

8 

0 

82 

79 

71 

58 

41 

21 
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31 

32 

29 

23 

16 

9 

0 

83 

80 

72 

59 

42 

21 

0 

34 

33 

29 

24 

17 

9 

0 

84 

81 

73 

59 

42 

22 
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35 

34 

30 

25 

18 

9 

0 

85 
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60 

42 

22 
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35 

31 

25 
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9 
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86 
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43 

22 
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84 

75 
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44 

23 
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38 

37 

33 

27 

19 

10 

0 

88 

85 

76 

62 

44 

23 
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39 

38 

34 

28 

20 

10 
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89 

86 

77 

63 

44 

23 

0 

40 

39 

35 

28 

20 

10 
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90 

87 

78 

64 

45 

23 

0 

41 

40 

36 

29 

20 

11 

0 

91 

88 

79 

64 

46 

24 

0 

42 

41 

36 

30 

21 

11 

0 

92 

89 

80 

65 

46 

24 

0 

43 

42 

37 

30 

22 

11 

0 

93 

90 

81 

66 

46 

24 

0 

44 

42 

38 

31 

22 

1 1 

0 

94 

91 

81 

66 

47 

24 
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45 

43 
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22 
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95 
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12 
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96 
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83 

68 

48 

25 
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24 
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49 
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42 
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24 
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35 

25 

13 

0 
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87 
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50 

26 

0 
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of  F(C)^,/F(C)^,  will  not  differ  appreciably  from  F(C)q/F(G)q.  Hence,  if  a 
scale  is  assigned  to  the  observed  F(hkO)  values  so  that  for  the  trial  structure 

used,  the  ratio  J |^0bs|  locale  | is  unity  *or  t^le  8rouP  CEFG,  this 

ratio  will  decrease  for  each  successive  group  like  CEFG,  and  from  the  numeri- 
cal values  of  the  ratios  a temperature  factor  can  be  calculated  (Section  11)  by  a 
simple  graphical  method. 

9.  Graphical  Refinement  of  a Trial  Structure 
The  structure  factors  in  the  group  CEFG  will  be  relatively  insensitive 
to  the  exact  location  of  the  atoms,  and  the  sensitivity  will  increase  with  each 
succeeding  group.  If,  therefore,  an  approximate  structure  is  postulated  to  give 
reasonable  agreement  with  the  F(hkO)  values  in  the  first  group,  it  can  be  im- 
proved by  comparison  with  each  group  in  turn,  in  the  following  way:  Suppose 
that  for  a projection  with  an  atom  located  so  that  its  (x,  y)  circle  is  centered 
at  Qj,  the  F(hk)  values  are  calculated  and  compared  with  the  observed  values 
(Fig.  11).  Suppose  that  agreement  would  be  improved  if  the  values  of  Aj,  B^, 

Cj.  . . were  altered  as  shown  by  the  arrows  to  positions  A^,  B^,  • • • The  best 
circle  that  can  be  drawn  through  A^,  B^,  C . . is  centered  at  Q^.  which  then 
represents  the  improved  position  for  this  atom. 


Fig.  11.  Process  of  refining  a projection. 
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10.  General  Utility  of  the  Construction 
The  construction  appears  to  have  some  value  in  the  fact  that  it  offers  a 
very  easy  way  of  visualizing  the  contribution  from  an  atom  oi  group  of  atoms 
to  a large  number  of  structure  factors  simultaneously  without  too  much  over- 
crowding of  the  diagram.  Space-group  extinctions  can  be  readily  represented, 

e.g.,  in  Fig.  12.  Consider  the  face-cen- 
tered [oo l] projection  in  which  the  atom 
at  the  origin  O is  represented  by  the 
circle  through  A,  and  the  atom  at  the 
face  center  by  the  circle  through  B. 

These  circles  cut  the  F(hk)  line  at  P and 
Q respectively,  and  OQ  = OP/2. 

In  Section  3 it  was  shown  that  OP  has 
to  be  graduated  so  as  to  have  (h  + k)  equal 
parts  with  maxima  at  O and  at  P.  Hence, 
if  h + k = 2n,  Q will  be  on  a maximum,  while 
if  h + k = 2n  + 1 , Q will  be  on  a minimum,  and 
atom  B will  cancel  out  atom  A's  contribution.  Hence,  the  space  group  condition  for 
a centered  [00l]  projection  can  be  derived;  namely,  reflections  are  present  only 
for  h + k = 2n. 

ii  kt  , „ , , /,  (-2B  sin2e)A? 

11.  A Nomogram  for  Evaluating  f^/f^  = e 

The  nomogram  has  4 scales  A,  B,  C,  and  D.  A and  D art.  logarithmic 
scales  with  1 cycle  = 10  inches;  they  are  equidistant  from  C,  which  is  the  log- 
arithmic scale  having  2 cycles  = 10  inches  required  for  multiplication. 

Since 

_ - ( -2B  sin^  0)/k2 

fT  = foe 

2B  sin^O 

-ng^rr  72— 

= f 10  e K 

o 


Fig.  12.  Demonstration  of  space 
group  extinction  conditions  for  a 
face-centered  projection. 


Fig.  13.  Nomogram  fcr  evaluating  f_ /f 
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the  D scale  has  its  origin  displaced  to  the  left,  so  that  the  reference  line  carrying 

the  origins  of  the  A,  B and  C scales  cuts  the  D scale  at  the  value  1.15.  Thus,  if 

the  value  of  (sin  0)/  X is  found  on  the  A scale,  and  the  value  ol  B on  the  D scale, 

r ^ ? 1 

the  line  joining  these  points  will  cut  the  C scale  at  the  point  I (B/ 1 . 1 o)  x (sin‘  '0)A  j- 
The  antilogaritbm  of  this  value  is  now  required,  and  the  B scale,  which  is  a 
linear  inch  scale,  performs  this  operation  against  the  A scale.  Note  that  the 
values  on  the  B scale  increase  from  right  to  left,  because  the  original  equation 
contains  the  exponent  with  a negative  sign.  An  example  of  the  use  of  the  nomo- 
gram for  evaluating  f^,/ f for  particular  values  of  B and  (sin^  0)/X^  has  h< .an 
sketched  in  on  the  diagram  (Fig.  13):  if  sin^ 0/X^  = 0.04  and  B = 3,  f T/fQ  = 0.79. 

This  nomogram  can  also  be  used  as  follows  to  evaluate  the  temperature 
factor  for  the  crystal  in  conjunction  with  the  procedure  described  in  Section 
When  a series  of  values  of  the  ratio  2 | Fobs K 2 |Fcaic|has  been  s.lo v-ted 

2 2 r 2 21 

for  groups  having  a mean  value  of  (sin  0)/X  = [(sin  0)/X  the  value  of  each 

ratio,  when  found  on  the  A scale,  determines  a point  on  the  B scale  and  hence  a 

r 2 2 1 

point  on  the  C scale.  If  the  points  on  the  A scale  having  the  values  ^(sin  0)/X 
are  joined  to  ‘■heir  corresponding  points  on  the  C scale,  these  lines  when  pro- 
duced, should  converge  on  a particular  value  of  the  temperature  factor  B,  which 
is  thereby  determined. 

I am  grateful  to  Professor  A.  vonHippel  for  his  interest  in  this  problem. 


